Phase transitions in ferro-antiferromagnetic bilayers with a stepped 

interface 

D.P. Landau a '*, Shan-Ho Tsai a ' c , Thomas C. Schulthess b 

a Center for Simulational Physics, University of Georgia, Athens, GA 30602 
b Center for Computational Sciences and Computer Science & Mathematics Division, Oak Ridge National Laboratory, 

Oak Ridge, TN 37831 
c Enterprise IT Services, University of Georgia, Athens, GA 30602 



Abstract 

We have studied magnetic ordering in ferro/antiferromagnetic (F/AF) bilayers using Monte Carlo simulations 
of classical Heisenberg spins. For both flat and stepped interfaces we observed order in the AF above the Neel 
temperature, with the AF spins aligning collinearly with the F moments. In the case of the stepped interface there 
is a transition from collinear to perpendicular alignment of the F and AF spins at a lower temperature. 
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Magnetic properties of ferro/antiferromagnetic 
(F/AF) bilayers can be very different from those of 
free F and AF films. These coupled systems exhibit 
a shift in the hysteresis loop (exchange bias) and 
larger coercivity than the free F film [1]. Other 
effects of the F/AF coupling include the order in 
the AF observed above the Neel temperature [2] 
and the perpendicular orientation of the F with 
the AF moments[3,4]. Understanding these effects 
and the nature of the F/AF interfaces remains a 
challenge. 

In this paper we use Monte Carlo simulations 
to study the effect of the interfacial exchange and 
roughness on the magnetic ordering in F/AF bilay- 
ers. The work is motivated by recent experiments 
for Fe304/CoO multilayers [2] . We consider a ferro- 
magnetic (F) film coupled to an antiferromagnetic 
(AF) film where the lattice is coherent across the 
F/AF interface. Each film is a bec lattice, with lin- 
ear sizes L x — L y — L < 96 and 12 staggered (be- 
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cause of the bee structure) layers of classical spins 
I S r I = 1 , which interact via the Hamiltonian 
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where (r, r') denotes nearest-neighbor pairs of 
spins coupled with exchange interactions Jp = 
5J > in the F film, J A = -J < in the AF 
film, and Ji = — J at the F/AF interface. Spins in 
the AF film have a uniaxial single-site anisotropy 
Ka — J, whose easy axis is along the y axis. The 
demagnetizing field on the F film is modeled with 
a hard-axis (Kp = —0.5 J) along the z direction, 
which is perpendicular to the F/AF interfacial 
plane. No external magnetic field is applied. We 
use periodic boundary conditions along the x and 
y directions and free boundary conditions along 
the z direction. We model flat interfaces that 
are fully uncompensated as well as uniformily 
stepped ones that are compensated on average. 
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The stepped interface has 6, L, and one spin per 
terrace in the x, y, and z directions, respectively. 

The F and AF order parameters are the uni- 
form (m) and staggered (m s ) magnetization per 
spin, respectively [6]. m x and m y denote the two 
components of m in the interfacial plane. We per- 
form Monte Carlo simulations with Metropolis al- 
gorithm at fixed temperature T. Typically we dis- 
card 3 x 10 5 Monte Carlo Steps/site (MCS) for 
thcrmalization and then use about 2 x 10 5 MCS 
for averages. Whenever not shown, error bars in 
the figures are smaller than the symbol sizes. 

The T-dependence of m shown in Figs. 1(a) and 
1(b) plus finite-size analysis (not shown) indicate 
an F ordering transition at T c « 9.3J/ks- A free 
12-atomic-layer AF film undergoes a phase transi- 
tion at the Neel temperature [5] Tm ~ 2.2J/fcs. For 
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Fig. 1. Temperature dependence of m, m s , m x , and m y (see 
text) for L = 60, with (a) stepped and (b) flat interface. 

both flat and stepped interfaces, the spins on the F 
film orient predominantly in a direction collinear 
with the easy axis of the AF at high T (see Fig.l), 
even above Tm- This is an indication that there is 
still order in the AF above Tjv due to the coupling 
to the ferromagnct. As T is lowered in the case of 
the stepped interface the F spins switch to orient in 
a direction that is perpendicular to the AF spins. 
Our results suggest that the onset of this perpen- 
dicular orientation is very sharp and it occurs at 
a temperature below Tm- The z-components of m 
and m s are very small for all T. In the absence of 
the AF film, spins on the F film have global rota- 
tion symmetry in the x-y plane, which is their easy 
plane. The preferential orientations of the F spins 



observed either below or above Tm result from the 
exchange coupling to the AF film. 

Fig. 2 shows m s versus T for different L for 
stepped and flat interfaces. In the former case 
the decay to zero of m s at T « 2.2 J/ks becomes 
sharper for larger L (sec the inset), suggesting an 
AF phase transition at Tm- In contrast, with a flat 
(uncompensated) interface there is no finite-size 
dependency of m s near Tm, suggesting the absence 
of an AF phase transition there. 
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Fig. 2. Staggered magnetization as a function of tempera- 
ture for stepped and fiat interfaces. 
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